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Abstract 

One of the prime tools to search for new light bosons interacting very weakly with 
photons - prominent examples are axions, axion-like particles and extra "hidden" 
U(l) gauge bosons - are light-shining-through-a-wall (LSW) experiments. With the 
current generation of these experiments finishing data taking it is time to plan for the 
next and search for an optimal setup. The main challenges are clear: on the one hand 
we want to improve the sensitivity towards smaller couplings, on the other hand we 
also want to increase the mass range to which the experiments are sensitive. Our main 
example are axion(-like particle)s but we also discuss implications for other WISPs 
(weakly interacting slim particles) such as hidden U(l) gauge bosons. To improve 
the sensitivity for axions towards smaller couplings one can use multiple magnets to 
increase the length of the interaction region. However, naively the price to pay is that 
the mass range is limited to smaller masses. We discuss how one can optimize the 
arrangement of magnets (both in field direction as well as allowing for possible gaps 
in between) to ameliorate this problem. Moreover, future experiments will include 
resonant, high quality optical cavities in both the production and the regeneration 
region. To achieve the necessary high quality of the cavities we need to avoid too 
high diffraction losses. This leads to minimum requirements on the diameter of the 
laser beam and therefore on the aperture of the cavity. We investigate what can be 
achieved with currently available magnets and desirable features for future ones. 



1 Introduction 

Many extensions of the Standard Model (SM) predict a "hidden sector" of particles which 
transform trivially under the SM gauge group and therefore interact only very weakly 
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with the "visible" sector, i.e. the particles of the SM. Due to their feeble interactions 
such particles are notoriously difficult to detect and even very light ones, with masses 
below an eV, may have escaped detection so far. Such new very weakly interacting sub- 
eV particles have been dubbed WISPs (weakly interacting slim particles). Beyond this 
purely phenomenological reasoning we have good theoretical motivation for the existence 
of light but also very weakly interacting particles. The prime example is the axion which 
arises as a very light pseudo-Goldstone boson in the course of a solution of the strong 
CP problem [IHZ]- Here, because of the pseudo-Goldstone nature, the weakness of the 
couplings ~ 1/ fa and the smallness of the masa3 rn a ~ m^f^j ' f a ~ meV(10 10 GeV// a ) are 
inherently related to a high energy scale f a at which the breaking of the Peccei-Quinn 
symmetry occurs. So, not only it is very plausible to have light, very weakly coupled 
particles, but indeed if we find them we may obtain information on underlying physics at 
very high energy scales. 

Similarly other light (pseudo-) scalar bosons with analogous couplings but possibly dif- 
ferent masses, so called axion-like particles (ALPs), arise quite naturally. Indeed top-down 
models such as those arising from compactifications in string theory suggest plenty of 
ALPs [8HTT] . Other particularly well motivated WISPs candidates are hidden U(l) gauge 
bosons [12l[T3] and mini-charged particles [13J. They naturally arise in models based on 
supergravity or superstrings [MHPQ] . All in all, the detection (or also the non-detection) 
of WISPs may give us important information on the underlying structure of fundamental 
physics. Although we keep an eye on general WISPs, the main focus of this paper will be 
the optimization of experiments searching for ALPs. 

Optical precision experiments are a powerful tool to search for WISPs, thereby exploring 
the hidden sector. In particular, laser experiments have an enormous potential to search 
for particles with tiny couplings to photons. A class of very simple and effective laser 
experiments is based on photon - WISP - photon oscillations: the so called light shining 
through a wall (LSW) experiments. Let us briefly recap the basic principle as shown in 
Fig. 1 (a) for the case of ALPs (Figs. l(b)|l(c) show the relevant processes for other WISPs) 



If light is shone through a magnetic region, a small fraction of the beam can be converted 
into ALPs, exploiting the Primakoff effect [2UH2"2] . A barrier after the magnetic regions 
allows only the ALP component of the beam to enter into a second identical magnetic 
region, where the ALPs can be reconverted into photons. 

The sensitivity of these experiments has grown considerably over the last few years, 
to the point that by now they are the most sensitive purely laboratory prob e£0 for ALPs 
in the sub-eV mass range. This is shown in Fig. [2] which displays the present status of 
constraints on the two-photon coupling g of ALPs versus their mass m^. LSW experiments 
currently reach a sensitivity g ~ few x 10 _7 GeV~ 1 . Clearly, the goal (see, e.g., J32j[33] ) of 
the next generation of LSW experiments should be to reach and maybe even supersede the 
present limits arising from astrophysics and from the non-observation of photon regener- 



1 m 7r , fir are the known pion mass and decay constant. 

2 Indeed laboratory bounds are less model dependent. They also apply if the couplings to photons 
effectively depend on environmental conditions such as the temperature and matter density 
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Figure 1: Basic principle of an LSW experiment. An incoming photon 7 is converted into a new particle 
which interacts only very weakly with the opaque wall. It passes through the wall and is subsequently 
reconverted into an ordinary photon which can be detected. The different panels show the explicit processes 
contributing to LSW for various WISPs. From left to right we have photon - ALP, photon - hidden photon 
and photon - hidden photon oscillations facilitated by MCPs. 

ation of solar ALPs in the magnetic field of the CERN Axion Solar Telescope CAST (cf. 
Fig. [2]) [M]. This requires a significant improvement in sensitivity by four orders of mag- 
nitude to the ~ 10~ n GeV -1 level. Beyond entering untested parameter space, ALPs with 
such a coupling are particularly well motivated, both from theory as well as phenomenol- 
ogy. Firstly, massless ALPs, with coupling to photons in the g ~ a/M s ~ 10~ n GeV -1 
range could occur naturally in string compactifications with an intermediate string scale 
M s ~ 10 9 GeV. Secondly, there are a number of puzzling astrophysical observations which 
may be explained by the existence of ALPs with g in the above range [35tf4"T] (the inter- 
esting areas are marked orange in Fig. [21 ) 

The most straightforward way to increase the sensitivity for ALPs is to enlarge, BL, 
where B stands for the magnitude of the magnetic field and L its length. This directly 
increases the probability of 7 — > ALP — > 7 conversion, which at small masses scales 
as oc {BLf. Nearly all the current LSW experiments (ALPS [SUET], BFRT [52l[53], 
BMV [5HE5], GammeV [56], LIPSS [SZ1ISH] and OSQAR [59]) performed so far, recycle 
one or two of the long superconducting dipole magnets from accelerator rings, such as the 
ones from HERA, Tevatron or LHC. A straightforward improvement can be achieved by 
using a larger number of these magnets (e.g. there are about 400 superconducting dipole 
magnets in the decommissioned HERA ring |60j). 

Another important step in advancing LSW experiments is the introduction of matched 
optical resonators in both, production and regeneration regions [HHE2], thereby increas- 
ing the probability for a photon to arrive at the detector by a factor ~ N paBB! xN pasSt 2, 
corresponding to the number of passes N pasSti in each cavity. 

Implementing these advances, sensitivities in the g ~ 10~ n GeV -1 range, for light 
ALPs, seem achievable, thus opening great opportunities for discoveries. However, both 
the increase in the length of the interaction region as well as the resonant cavities come at 
a price. An increased length of the interaction region actually reduces the ALPs sensitivity 
at larger masses, and high quality cavities must have a minimal diameter in order to avoid 
excessive diffractive losses, in turn requiring magnets with sufficiently large aperture. It 
is the purpose of the present paper to discuss these effects and ways how to ameliorate 
them. Moreover, we want to find an optimal configuration, based on existing technology, 
that maximizes the sensitivity of the experiment for a wide range of ALPs masses. 
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Figure 2: Summary of cosmological and astrophysical constraints for axion-like- particles (two photon 
coupling g vs. mass m^, of the ALP). Areas with interesting astrophysical hints, e.g. the one for a non- 
standard energy loss in white dwarfs |43| or the one for an anomalous 7-ray transparency of the universe 
(e.g. [44[[45]), are marked in orange. The parameter range for the axion is shown hatched. Note that the 
limit from the microwave cavity axion dark matter search experiment ADMX [48] is valid only under the 
assumption that the local density of ALPs at earth is given by the dark matter density. (Compilation from 
Ref. [49], where also details can be found.) 

The paper is set up as follows. In the following section we briefly recapitulate the 
necessary formulae for LSW experiments. In particular we review the basics of photon - 
WISP - photon oscillations and the effects of resonant cavities. In Section [3] we discuss how 
the loss of sensitivity at large ALP masses can be ameliorated by alternating the direction of 
the magnetic field (see also Refs. [22lf63] ) and leaving suitable gaps in between the magnets 
- a previously unnoticed, but crucial feature which can be exploited, too. Section H] is 
devoted to the diffractive losses in the optical cavity and the resulting requirements for 
the aperture of the magnets (see also [61] ) . We then study what is possible with currently 
available magnets and what would be ideal features of future magnets. Finally, in Section [5] 
we summarize our findings and conclude. 
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2 Basics of photon- ALPs-photon oscillations 



In this section we would like to recapitulate the theoretical basis of LSW experiments, 
in particular to recall the general form of the oscillation probability and then deduce its 
explicit expression for the setup of current experiments. 

At low energies, the coupling of an ALP, 0, to photons can be described by the effective 
Lagrangian 

C = - l -F^F, u + 1 {d^dy - m|0 2 ) - \g^ v F^, (1) 

where F^ u is the electromagnetic tensor and F^ v = i>€ l j lV p\F pX its dual. This leads to the 
following equations of motion [65], 

(2) 



where n stands for the refractive index of the medium through which the photon propagates. 
At first order in gBL, with L the linear dimension associated with the extent of the 
magnetic field [2U1I22], one finds two possible solutions, 

^V- 6 gfi(f )■§(?), (3) 

Air r — r 



where = y u 2 — m^. Specializing to the experimentally relevant case that the photon 

beam is sent along the x-axis and that the transverse extent of the magnetic field is much 
larger than that of the laser, the problem becomes one- dimensional [22] . 

<j>(x,t) = e -i^- k ^)^- f dx'Eix 1 ) ■ B(x'). (4) 

J 

This can be further evaluated by inserting the appropriate plane wave form E(x,t) = 
e 2 i?oe iaJ ^ na;_ ^ for the electric field of the laser beam, assumed to be linearly polarized in the 
^-direction, and expressing the magnetic field as B(x) = e z B f(x), where f(x) accounts 
for the spatial variation of the field, 



[x 



t) = ^-Bue'^'^) J dx'f(x') e iqx ', (5) 



where 



m 2 



q = nu — \ ui 2 — m 2 m oo(n — 1) + — - (6) 

is the momentum transfer to the magnetic field. Therefore, the probability that a laser 
photon converts into an axion (and vice versa) after traveling a distance L can be written 

as 

P^ = P^ = (gBL) 2 \F(qL)\ 2 , (7) 
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where the function F(qL) is the so called form factor, defined as 



F(qL) 



1 

I 



dx'f(x') 



D iqx 



So far for the general description of laser photon <H- ALPs oscillations along an arbitrary 
transverse magnetic field. Now we want to specialize to the case exploited in the current 
generation of LSW experiments, exploiting single dipole magnets with a homogeneous field 
along the z-axis, both on the generation as well as on the regeneration side. In this case, 
the explicit expression of the form factor is 



| -^single (<?£)| 



2 . fqL 
— - sin — 
qL V 2 



(9) 



Here, the subindex "single" is meant to stress that this expression is valid only for a single 
homogeneous magnetic region. In this case, the maximum conversion probability, 



is achieved in vacuum 
small masses, 



n 



« g 2 B 2 L 2 /A, (10) 
1) for small momentum transfer, qL/2 <C 1, corresponding to 



< 8.9 x 10~ 4 eV 





1/2 


"nr 


1/2 


LeV. 




X 





(11] 



where the form factor takes its maximum value, F s - mg \ e (qL) ps 1. Correspondingly, the 
constraints on ALPs from LSW experiments are best at small masses, cf. Fig. [2j For 
larger masses, the oscillations of Eq. (jUJ) start to dominate, and the (envelope of the) upper 
limit on the coupling starts to grow as m^, represented in the double logarithmic plot Fig. [2] 
as an approximate straight lin^E Given the formula (jZJ) it is clear that efforts should be 
made in the setup of LSW experiments in order to keep F(qL) ~ 1 as much as possible, 
for a wide range in <C u. 

Finally, we would like to note that LSW experiments are sensitive not only to ALPs, 
but also to other WISPs, notably to hidden-sector U(l)'s and mini-charged particles 



(MCPs) (cf. Figs. [T(b)1 and [T(c)j). 

The interactions of the former with the standard photon are described by the low energy 
effective Lagrangian 



C D 
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1 fu/ 1 



1 Y 

-X X^ v - —X F^ v 



^m 2 Y X,X», 



(12) 



where A^X^ are the electromagnetic and hidden vector potentials, F^ v = d^A y — d v A 
and X^ = d^Xy — d u X^ their respective field strength, x 1 is the kinetic mixing 
parameter [13], and my is the hidden photon mass. Similar to the ALPs case, kinetic 

3 For axion-like particles with masses of the order of the frequency of the laser ~ uj, a resonant 
oscillation appears and the above formulas are not valid anymore [66] . 
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mixing gives rise to photon - hidden photon oscillations. The corresponding probability 
can be obtained by solving the respective equation of motion, 



(u, 2 + dl) 1 - 



-2u 2 {n- 1) X 2 

-x 



-x 



0. 



and reads 



P. 



7— >7' 



4x sin 



qh 



tot 



where q is now given by 



nuo 



uj 2 - m 2 , 



ui(n 



in 



7' 



2u 



(13) 



(14) 



(15) 



and where L tot is the total distance of travel. Clearly, for photon - hidden photon oscilla- 
tions no magnetic field is needed. 

Finally, the relevant equation of motion for the case of photon - massless hidden pho- 
ton (the latter kinetically mixing with the photon) oscillations via a loop of minicharged 
particles (cf . Fig. 1 (c) ) reads 

{u 2 + d 2 z )i ' «■ 



2cu 2 e 2 h AN l 



-x 



0. 



(16) 



where is the unit hidden sector U(l) charge, % =||,-L indicate the polarization with 
respect to the magnetic field, and AiVj are the magnetic field dependent, complex refractive 
indices describing the refraction and absorption due to the virtual and real production of 
MCPs. See Ref. [67] for details, where also the explicit expression for the corresponding 
oscillation probability can be found. 



3 Refinements of LSW experiments 

For the next generation of LSW experiments several important proposals have already been 
made to improve the sensitivity. One of the most effective novelties is to include matched 
Fabry-Perot cavities in the production and regeneration sides of the experiment [6TII62] . 
When both cavities are tuned to the same frequency, co, it is possible to gain an enhance- 
ment in the sensitivity for the ALP-photon coupling or the kinetic mixing parameter by 
the fourth root of each of the cavities' power buildup^, g scns or x scns oc {f3 g f3 r )~ l / 4: . Con- 
sidering that with the available technology cavities with (3 ~ 10 4 — 10 5 seem realistic, an 
improvement of the order of 10 2 in these couplings are feasible. The expected number of 
photons after the regeneration cavity of such an LSW experiment will be given by [6T|[70] 

N s = V 2 (3 g (3 r 3™ P^^t, (17) 

4 From now on, we chose to refer to the power buildup of the cavity, instead of the commonly used 
finesse since it is the former which plays the most direct role in the production and regeneration of WISPs. 
See also Sec. |U 
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6 + 6 HERA magnets U = 8.8 m) 


L = 52.8 m 


Magnetic field 


B = 5.5 T 


Primary laser power 


V prim = 3 W 


Power build-up 


P g = /3 r = 10 5 


Laser frequency 


u = 1.17 eV 


Overlap between WISP mode and electric field mode 


77 = 0.95 


Detection time 


r = 100 h 


Dark count rate 


n b = 10" 4 Hz 



Table 1: Benchmark values for a next generation LSW experiment. 



where (in a symmetric setup) P 7 ^wisp — -Pwisp^7 is the probability of photon-ALP con- 
version, e.g. Eq. ([7]) or (JHJ), P pr im is the primary laser power, /3 9 r are the power build-ups 
of the generation and regeneration cavities, rj is the spatial overlap integral between the 
WISP mode and the electric field mode [61] and r is the measurement time. 
The significance of the WISP discovery scales as S 



2(^N S + N b - v^Vfe) where N b 
is the expected number of background events [7T|[72]. To estimate the sensitivity of a 
LSW experiment let us assume that no significant signal over background is found. Then a 
95% CL. exclusion limit can be set requiring S < 2 which translates into N s < 1 + 2\fW b . 
Writing N s = n s T and N b = n b r with n s and n b the detection rates of signal and background 
photons we find that the limit on n s scales improves with measurement times as 1/r 
if backgrounds are negligible and only as 2^n b /r where the experiment is background 
dominated. Let us consider that we are in the second situation^]. The corresponding 
expected sensitivity for the ALP-photon coupling of such an experiment is thus 



2.71 x 10 



-ii 



GeV 



\F(qL)\ 



290 Tm 



0.95 
L V 



1/2 



"10 10 " 


1/4 


"3 W" 


1/4 


■ n b ' 


1/8 


"I00h" 










ao- 4 Hz. 




T 



where we have used the benchmark values for the most important parameters, as summa- 
rized in Table [TJ 

Since an improvement in the length of the magnetic region seems mandatory for future 
LSW experiments, we have assumed a benchmark scenario with 6 superconducting dipole 
magnets on the production and the regeneration side of the experiment (cf. Tab. [T]). 
Correspondingly, the longitudinal profile f(x) of the magnetic field will not be so simple as 
in the single magnet setups of the first generation of LSW experiments. Therefore, the form 
factor F(qL) in Eq. ( Tl8|) has to be determined via Eq. ([8]) using the actual longitudinal 
profile. It is the purpose of the remainder of this section to do this for realistic arrangements 
of magnets and to determine its actual influence on the sensitivitjj^]. 



5 If feasible, it makes sense to extend the measurement period until the 1/y/r regime is reached. 

6 Clearly, the sensitivity to the kinetic mixing between the photon and a massive hidden photon will 
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3.1 Generalization to a series of magnets with gaps in-between 



The setups foreseen for the next generation of LSW experiments will exploit series of N 
dipole magnets, including a natural and probably unavoidable "gap" with no magnetic field 
between each magnet. Therefore, let us investigate the form factor (jSJ) for a longitudinal 
profile corresponding to iV equally spaced magnets, each of length £, separated from each 
other by a fixed length A. In fact, a short calculation results in 



with L = N£ the total length of the magnetic field. 

The interpretation is clear. The form factor is the sum of the individual form factors 
of the production in the different magnets (which being equal factorizes out), weighted by 
a phase that accounts for the phase delay of the ALPs produced in the different magnets 
(Since a common phase is arbitrary, we have measured these phases with respect to the 
axions produced in the first magnet, so it is e ii ? n ( £ + A ) for the (n + l)th magnet). 

This function has two types of zeros. First, it is evident the zeros of the form factor of 
one single magnet (q£ = 2nk, k £ Z + ) should be still there since it appears as a common 
factor. Secondly, there are the zeros related to the gaps, appearing when qN(£ + A) = 
2nk', k! £ Z + and q(£ + A) ^ 27r/c", k" £ Z + . For non-zero A these conditions are fulfilled 
when k' is not an integer multiple of N, i.e. k' = 1, 2, N —1,N + 1, ...27V - 1, 27V + 1, 
etc. This means that for each zero of the single magnet form factor there are N — 1 zeros 
coming from the gaps. When A — > the combination of the two types of zeros coincide 
with the zeros of a magnet of length L. As a check, we note that Eq. (1191) approaches the 
previous single magnet result (J9]) in the limit A — > or N = 1, as it should. 

All in all, the presence of gaps does not alter the number of zeros of the form factor, 
which ultimately means dips in sensitivity we want to avoid. It does however modify the 
values of q and hence the ALP mass where they appear. As we increase A from the zeros 
displace towards smaller values of q (and hence of m$). In particular, the first zero of the 
gapless form factor in Eq. (Q, q = 2n/L, now moves to q — 2ir/(L + TV A). This means 
that ultimately, the price to pay for enlarging the magnetic region by piling up individual 
magnets with gaps in between is the loss of the coherence of the photon-ALP conversion 
already at smaller masses (with respect to the ideal gapless case). 

Not surprisingly, our result coincides with that obtained by introducing phase shift 
plates inside the magnetic region [73]. In [73], the authors argue that it is possible to put 
the photon and ALP waves in phase after some distance by the insertion of a phase shift 

not depend on the longitudinal profile of the magnetic field. Still, as we will see, the laser improvements 
and the increased length of the generation and regeneration sides will lead, in the next generation of LSW 
experiments, to a sizable improvement in sensitivity. 




e iql _ 1 \ 1 _ e iqN(i+A) 

qL ) 1 - e «z(<+A) 



(19) 
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0.53 



0.55 



2.51 
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0.92 



Figure 3: Sensitivity for an LSW experiment using 6 + 6 HERA magnets. The considered benchmark 
parameters of the setup are given in Tab.[TJ The enhanced zone corresponds to = 4.7 — 5.5 x 10 -4 eV, 
and the size of the respectively needed gaps are shown in the table. The fc's are the corresponding integers 
in Eq. (|20l) . 



plate that advances or delays the photon phase with respect to the axion wave. In our 
case, we can do the same after each magnet by choosing the appropriate gap length^. The 
clear advantage is that the gaps do not introduce extra optical losses. 

Since the gaps' length appears as a free parameter which is adjustable to our needs 
(besides the number of magnets, which can not be altered significantly), we can optimize 
the form factor Eq. ffl9|) for a given non-zero ALP mass by choosing the appropriate value 
of A. The maximum of Eq. f fT9|) happens when the ALP production in all the magnets 
adds up in phase, that is when the phases in the sum in Eq. ffl9|) add up coherently, i.e. 
when 

|( 1 + j)=^ ( fceZ+ )- ( 2 °) 

At these maxima, the form factor is reduced to the one corresponding to a single magnet 
of length £, i^AUa* = F(q£)- This is the best one can achieve with a multimagnet 
configuration without using phase-shift plates or filling the oscillation region with a buffer 
gas [51]- 

However, the fact that this optimization is mathematically possible does not mean 
that it can be conveniently implemented in a realistic experiment. As we can easily see, 
if we want to optimize for instance the low mass region q£/2 < 1 we shall have at least 
q£(l + A/£)/2 = 7r, requiring A > £, which is a very large gap. In Sec. |4]we will study the 
limitations on the length L imposed by the current available magnets and will understand 
that the total length is a valuable goocfl 



7 We note, however, that this is only true for ALPs. For hidden photons, which do not interact with 
the magnetic field, it is only the total length (including the gaps) that counts. 

8 Indeed in some cases it might even be useful/possible to fill the unavoidable gap with another magnet. 
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Table 2: First maxima of the form factor for an array of magnets alternating in n subgroups. 
For n — > oo, qL/2 — > rar/2. 



n = 


1 


2 


3 


4 


5 


6 


qL/2 = 





2.331 


4.131 


5.832 


7.486 


9.115 



We conclude that changing the gap size allows for filling the sensitivity dips of the 
form factor, improving the sensitivity of LSW experiments especially at large masses cor- 
responding to q£/2 > 1. For the benchmark setup proposed in Tab. [H this allows to 
significantly improve the sensitivity in the mass range ~4x 10 -4 — 10 -3 eV. If we try 
to extend this down to smaller masses, the size of the gap becomes excessively large for a 
realistic experimental setup. For instance, in order to maximize the sensitivity for an ALP 
mass of = 10~ 4 eV, an enormous gap between the HERA type (£ = 8.8 m) magnets 
of A = 285 m would be needed. As an illustration, we show in Fig. 13.11 the projected 
sensitivity for a 6 + 6 HERA magnet configuration, taking into account gaps between the 
magnets. The enhancement corresponds to an ALP mass range of (4.7 - 5.5) xlO -4 eV. 
The input details of the configuration are the same as the ones shown in Eq. ffl8|) . We 
have chosen values of A that are reasonable experimentally. 

3.2 Generalization to a series of magnets with alternating polar- 
ity 

Now we will extend this study considering an array on N identical magnets of length £, 
segmented into n subgroups of alternating polarity, such that the total magnetic length is 
given by L = N£. This idea was proposed in Ref. [22] as a technique to peak the form 
factor at nonzero values of q. The formula, neglecting the gap between the magnets is 
given by [22|l63] 

\F n (qL)\ = 

The maxima of this form factor (with respect to qL/2) happen close to the poles of the 
tangent, which are at the same time zeros of the accompanying sine or cosine. The pres- 
ence of the external factor 2 / (qL) slightly moves the maxima away from the poles of the 
tangent, this correction disappearing at large n (corresponding to large qL/2). The abso- 
lute maximum is the one with smallest qL/2. We have searched numerically the absolute 
maxima of the form factor in Tab. [2] for different values of n. For sufficiently large n we 
can treat the the problem analytically performing an expansion in l/(qL). One can see 
that the solutions approach asymptotically qL/2 = mr/2, and for them F n (qL) — > 2/tt. 



A S in (f) tan (g) 



qL 



cos (f ) tan (£) 



// even. 
n odd. 



(21) 



11 



If we now include a fixed gap in-between the magnets, as we did in Eq. (fT9|) . we find 
\F N:nA (qL)\ = 



n /gL\ sin(^(L/jV+A)) (qN_(L_ , an\ 

/gL\ cos(^(L/7V+A)) r?JV L , 



qL Sm \2N) sin(f(L/7V+A)) 

In the limit A — > 0, formulas (T2T]) are recovered. Moreover, the non- alternating form factor 
(Eq. (fT9|) ) is recovered for n = 1. 

Depending on the choice of the pair {N, n}, there is a different maximum for the form 
factor. Indeed, they can be searched as before by finding the poles of the tangents. In 
general we find that they are close to Nq(£ + A)/2 = n(2k + l)vr/2, k G Z + and the most 
important one corresponds to k — 0. This is asymptotically true as Nq(£ + A)/2 becomes 
large but again, the poles have to be searched for numerically for 0(1) values. In general 
and as with the previous case, the maxima displace towards smaller q when we increase A. 

We can again use the gap length A to optimize the form factor for a given ALP mass. 
The absolute maxima of Eq. ( 12"2"|) cannot be found analytically except in particular cases, 
however it is easy to see that they lie either very close to the poles of the tangent or to the 
zeros of the sine in the denominator (only for the n-odd case). The poles of the tangent 
are 

— (£ + A) = - + 7ik, keZ + , (23) 
2n 2 

and for these values the form factor evaluates to 

nlN 

F N , nA {qL) = F smglc (qi) - - y (24) 

S111 TV U + 7TK ) 

The zeros of the sine in the denominator are given by Eq. ( 1201) (they appear also in the 
non-alternating configuration) and for them the form factor evaluates to -Fz^n-odd.A^-k) = 
F s i ng \ e {q£)/n. 

The N = n case (i.e., a titi ••• configuration) can be solved analytically. One finds 
that the absolute maxima coincide exactly with the poles of the tangent. In this particular 
case, full coherence between the ALP production in all the magnets can be achieved and 
correspondingly the form factor saturates to the single magnet form factor F^^ A (qL) = 
F s i ag i e (q£) (In the second set of maxima, Fn,n,a is suppressed by a factor 1/N). 

The case n = 1 is trivial since it reduces to the non-alternating configuration. In this 
case, as we have said one can obtain also F N i A (qL) = F sing i e (q£) but this happens for the 
zeros of the sine in the denominator and the poles of the tangent produce non-absolute 
maxima for which the form factor is a bit smaller (is easy to see that Fn,i,a{qL) — > 
F s - mg ie{q£) x 2/V in the N — > oo limit). 

All in all we see that for certain ALP masses, the sensitivity can be nearly fully restored. 
Also, it is easy to see that this optimization can be performed even in the case qljl < 1, 
unlike the non-alternating configuration considered in the previous section. As an example, 
in the case of iV = n HERA magnets, with gaps of A = 1 m, the first maximum condition, 
k = 0, is reached for = 3.87 x 10~ 4 eV, a region where the previous setup (without 
alternation and gaps) had no sensitivity, cf. Fig. 13.11 



(22) 

n odd. 
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Figure 4: a) "6+6" configuration of HERA magnets. Combining all the possible wiggler configurations it 
is possible to restore the sensitivity of the experiment until the first zero of the form factor, where qL/2N = 
ir/2. The size of the gaps are A[m] = {4.5,3.9,3.4,2.5,1.9,0.9}. b) Sensitivity of a toy configuration of 
"400 + 400" LHC magnets in order to touch the predictions for the QCD axion ("axion band"). Using a 
gap in-between the magnets of A = 0.6 m we are able to scratch the axion band in ~ 3.15 x 10" 4 eV. 
In order to scan an ALP mass of ~ 2.4 x 10 -4 eV, a gap of A = 11.2 m is needed. However, in order 
to get the sensitivity shown in this plot the bore aperture radius of the magnet should be at least 16 cm. 



3.3 Some examples 

As an example, let us consider a possible 6 + 6 HERA configuration, i.e. 6 magnets in the 
production and in regeneration side, respectively, with the benchmark parameters given in 
Tab.H 

With six magnets on each side of the experiment we can produce four different symmet- 
ric configurations: tttttt> T4T4T4, TT44TT, TTT444- In Fig- H a) we display the sensitivity 
of all four possible 6 + 6 configurations, taking into account seven different gap sizes. We 
infer that, using these alternating field configurations, we are able to restore the sensitivity 
almost completely up to masses of order ~5 x 10 -4 eV. 

The gaps in the sensitivity, where qL/2N = kn, with k G Z + , can not be covered with 
this setup, since all the form factors are zero there. However, the ALPS Collaboration [51] 
was able to fill such regions introducing a gas inside their cavity, changing the refractive 
index in the generation and regeneration zone. Another possibility to cover these regions 
is to use phase shift plates, or a movable wall |56j . 

An important issue for the next generation of LSW experiments would be to know their 
potential to achieve sensitivity in the parameter range of the proper QCD axion. In Fig. H] 
b) we show a (way too optimistic) toy example in order to answer this question: We have 
chosen an alternating field configuration, with N = n. According to Eq. (1231) . we select the 
appropriate gap to peak in the mass range «3x 10~ 4 eV. We used 400+400 magnets, 
with a magnetic field of 9.5 Tesla and a length of i = 14.3 m each. We see that even this 
overly optimistic configuration can just scratch the QCD axion! 

Clearly, this setup is unrealistic: it needs more than 400 LHC magnets (£ = 14.3 m, 
B = 9.5 T) in each the production and the regeneration side of the experiment, enlarging 
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waist Zr 

Figure 5: Collimated waist region of a Gaussian beam. The Rayleigh range is denned as the distance 
which the beam travels from the waist before the beam area doubles. It is the characteristic distance 
before the beam begins to diverge significantly. 

the size of the cavity to ~ 5960 meters! As we will see in the next section this entails, 
that in order to assume the benchmark values of Table [TJ magnets with an aperture bore 
radius of 16 cm are needed. This is considerably more than the 3 cm of the existing HERA 
magnets. 

In short, with the current designs LS W experiments will enter an unexplored region in 
parameter space, but are still far from being able to test the QCD axiorj§. 

4 Expected sensitivity for next generation of LSW 
experiments 

In the preceding section we have discussed the possibility to change the length between 
the magnets as a technique to improve the sensitivity of the experiment for certain ALP 
masses. Now we would like to make the discussion more realistic as well as investigating the 
possibilities for the next generation of photon regeneration experiments, given the actual 
equipment. For this, the maximum length of the cavity is crucial, since on the one hand 
the probability of photon-ALP conversion for small masses scales as L 4 , and on the other 
hand, bigger cavities will give us more freedom to choose gaps between the magnets. 

4.1 Diffractive losses 

Concerning the cavity, we can estimate the approximate optimum length for the different 
superconducting dipoles available. The bore aperture of the magnets is a key factor in 

9 Coherence is kept up to higher masses in LSW experiments exploiting more energetic photons, e.g. 
from synchrotron radiation sources or X-ray free-electron lasers |74 [ I75| . However, currently the sensitivity 
is far away from the QCD axion because the average photon flux of current or planned facilities is many 
orders of magnitude smaller than the one reachable with optical lasers (for a pioneering experiment in this 
respect see Ref. [7^). Besides, an important disadvantage is the impossibility to set-up cavities. 
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order to determine the length of the cavity. 

Let us consider a collimated Gaussian beam, propagating in the z direction. The point at 
which the beam has a planar wave front with curvature R = oo, is known as the beam 
waist, and it is characterized by a particular spot size, wq. The power intensity of a beam 
falls off very rapidly with the distance, and at some arbitrary length z, the spot size can 
be described by w(z), 



w(z)=w \ll+^—) • (25) 

The distance zr is known as the Rayleigh range and defines the point at which the 
beam area doubles. It is given by zr = ttwq/X, where A is the wavelength of the laser 
in the medium. The sketch of the configuration is shown in Fig. [5j The radial intensity 
variation of a Gaussian beam with spot size u is given by 

J(r) = ^e- 2 ^ 2 , (26) 

where V is the total power of the optical beam at the position of the waist. The fraction 
of the power transmitted through a circular aperture of radius a is given by 

J l[r)rdr 

The peak amplification factor of the incident power laser inside an optical Fabry-Perot 
cavity - known as power buildup (/?) - occurs at the resonance condition, uj = kir/L, where 
k G Z, and L is the length of the cavity, and is given approximately by [77J 



45 1 

(s + <5i + s 2 y 



where Sq accounts for the round-trip losses in the cavity and 8^2 are the losses due to the 
transmissivity of the cavity mirror^]. In order to maximize (3 we have to keep 8q and 82 as 
small as possible. Then the optimum value for 8\ is 8\ = 80 + 82 corresponding to a so called 
impedance-matched cavity (0 = (5o + ^2) -1 )- Let us remark that for an impedance-matched 
cavity, a simple relation can be established between the power build-up and the cavity's 
finesse 

T = ix0. (29) 

For a cavity in good vacuum conditions, one of the dominant contributions to 8q arises 
from clipping the cavity mode. If the beam waist is located near one of the mirrors, then 
5q 1ip = e~ 2 °" 2 l w2<yZ ^ where Z is the distance of the waist to the cavity end mirror. If the waist 
lies in the center the clippling occurs near both mirrors and a factor of 2 arises. We shall 



10 Thc formal definition of the delta coefficients is = e~ Si , where Ri stands for the reflectivity of the 
i-th mirror. The coefficient 80 is related to the internal losses inside the cavity, therefore <5o = Aa$L where 
is the attenuation coefficient in the cavity. 
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only consider the former case, which is more advantageous since it optimizes the spatial 
matching of the generation and regeneration cavity modes [BT] . Other contributions to 5 
like non-Gaussian beam profiles, imperfections in the mirrors, etc., can be subsummed in 
an additional parameter 5$. The laser frequency will be also crucial to achieve maximum 
sensitivity: if it is too large, dispersion in the mirrors might increase the roundtrip losses. 
For this reason we have chosen to use 1064 nm wavelength light all throughout the pa- 
per. The experience gained in the gravitational wave community in building high finesse 
cavities containing extremely intense beams with this type of light makes this choice most 
practical [75] . 

Let us now describe one possible optimization procedure for the design of a LSW 
experiment. We start with a number of dipole magnets of field strength B and aperture a 
and ask for the parameters of the cavity that optimize our sensitivity for the ALP coupling 
g. Let us focus in the low mass range, <C 10 -4 eV, larger masses will be dealt with later 
on. Our final sensitivity to the coupling constant g, given by eq. f llSp . depends mainly on 
the combination 

i i / _ \ 1/2 

e ^ 0+ w +S * + S2 \ (30) 

where we have used the impedance matching condition and substitute in terms of wo. 
We have taken the generation and regeneration cavities to be equal, (3 g = (3 r = (3, to 
maximize the bandwidth overlap. 

The above function can be easily minimized. First it is clear, that the highest sensitivity 
is reached for the lowest possible value for 5$ + 62- Currently the best mirrors are in the 
6-2 ~ 10 -5 ballpark so we will keep this number for the rest of the paper. We will also assume 
that 5q's smaller or of the same order can be achieved. We have now three remaining free 
parameters, L, Z and Wq. Optimization with respect to Wq is trivial, in order to render 
the exponential as small as possible the argument has to be maximum and this happens 
at Wq = ZX/ir (i.e. at z R = Z). 

In order to minimize with respect to Z and L we first substitute the condition = Z 
in Eq. p]). 

1/2 x 1/2 

(e- *r + 6* + 5 2 ) . (31) 

Note that Z > L since the length of the cavity Z has to encompass the magnetized length 
L = N£ plus the length of the gaps, NA, plus maybe a field free zone near the mirrors of 
the cavity, called dz. There is a certain difficulty in defining properly a minimization since 
Z and L are related but Z can grow continuously and L only in multiples of the magnet 
length I. To start with let us ignore these difficulties and just consider an L ~ Z which 
can be varied continuously, i.e. iVA + dz is small (we want to have as long a magnetic 
region inside the cavity as we can). We find the minimum with respect to Z to be given 
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implicitly byo 



9 

7ra 



2ZA 



S* + 5 2 . (32) 



This equation can be solved numerically but an approximate solution is very easy to 
find. Ignoring the parenthesis for a while, the equation has the solution Z = —ira 2 / (A log(<5Q + 
62))- This solution makes the parenthesis — log Sq + 82 — 1 ~ 4.75 (taking Z = L). The 
real solution will be a bit smaller such that the exponential kills this additional factor. As 
an example, for 8$ + 62 = 10 -5 , the simple estimate gives 0.08697ra 2 /A while numerically 
we find the optimum to be 

7ra 2 / a \ 2 1064 nm . . 

Z opt = 0.0755— = 89.2 m , 33 

A V20 mm/ A 

which corresponds to ixa 2 / (Z opt X) ~ 13.25 and then to 8q 1P /(Sq + #2) = 0.177. In previous 
literature [70] a value 8q 1P /(8 + 8 2 ) ~ 1 was used, leading to a slightly less optimal setup 
(see the dashed line in Fig. [6] (right)). 

Let us summarize our optimization procedure. First we have optimized with respect 
to 62, the smallest value always results in an optimum sensitivity. Then, the impedance 
matching condition gives an optimum for 8\ = 8q + 62- Since Sq encompasses the clipping 
losses it is a function of the length of the cavity and has to be maximized together with 
the total magnetic length. Equation ( 132]) tells us precisely the optimum cavity length Z opt 
taking into account these relation and has to be solved numerically by specifying a value 
for Sq + 82 and a relation between L and Z (assumed L ~ Z). 

Now we can tackle the fact that L can only grow by multiple integers of £ + A. This can 
be numerically computed. In Fig. E] (left) we have plotted as a black solid line the quantity 
to minimize (Eq. fl3"U|) ) as function of Z using L = Z for an array of LHC magnets for which 
a = 28 mm. The black dots represent the value of Eq. (130!) as a function of the number 
of magnets, i.e. Z = L = N£. The minimum of the continuous line stands clearly around 
Z ~ 175 m (our minimization formula, Eq. (|33|) . gave Z opt = 174.78 m). The possible 
realization closer to it corresponds to an array of 12 LHC dipoles. Note however that the 
minimum is reasonably wide and there is some freedom to chose iV without spoiling the 
sensitivity too much. 

Let us now move on to discuss the high mass range. The setup described above maxi- 
mizes the sensitivity at — > (equivalently q — > 0) but will have some unavoidable dips 
in sensitivity corresponding to the zeros of the form factor F^^(qL) of Eq. (fT9~l) . So we 
place ourselves in the situation where we have already performed a run in the previous 
mode and we want to scan further parameter space by a small modification of the set up. 
The dips of FN t &(qL) correspond to the unavoidable dips of the single magnet form factor 
-Psingie(^), i-e. q£/2 = kn, k E Z + plus the dips related to the sum of the different magnet 
contributions, i.e. qN(£ + A)/2 = jfcV, k' G Z + . 



11 As it turns out this equation even holds if we take Z = N(£ + A) and optimize with respect to the 
number of magnets, N. 
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To get rid of these dips we have four different techniques at our disposal, namely a) 
introducing phase shift plates, b) a buffer gas in the oscillation region, c) alternating the 
magnet polarity and d) changing the gap size. The first two techniques can be applied to 
optimize the dips of the single magnet form factor, i.e. can be used to optimize F sing \ e (q£) . 
The remaining two can never avoid these zeros. An optimal technique to cure all these 
dips in -F S i n gie(^) with an extra measurement with buffer gas in the oscillation region was 
described and implemented in [51]. With phase shift plates a better result is in principle 
possible but the optical dispersion caused by them (we would need in principle one per 
magnet) will affect the power build-up of the cavity more than the scattering on the rarefied 
gas required by the buffer gas technique, so in principle the latter is preferred. 

We are now left with the zeros produced by the interference of the ALP waves in the 
different magnets. At low masses, now meaning qijl < it, i.e. below the first zero of the 
single magnet form factor we can improve substantially by performing our ALP search 
with the different configurations alternating the magnets. As we have seen, this procedure 
allows almost to maximize the form factor (one can get F ~ 2/ir ) around certain values 
of qijl in the low mass region, cf. for instance Fig. HJ 

We can finally play with the size of the gaps. As we have seen, they allow to shift the 
pattern of zeros and maxima of the form factor. However, we find ourselves between a rock 
and a hard place. On the one hand we would like to have as much magnetized region as 
possible (respecting the minimization of g SC ns), for instance in the above example we would 
chose A = 12 instead of N = 10 LHC dipole magnets. On the other hand, we would also 
like to have as much gap space to widen as much as possible the sensitivity peaks. A mixed 
solution would be to perform a first search with N = 12 and then change to N = 11 or 
N = 10 without changing the cavity length, such that the range of available gap sizes is 
bigger. 

In Fig. E] (right) we show how the different set ups loose their maximum sensitivity 
as we increase the gap size. When the total length reaches the 160 ~ 170 m figure, the 
sensitivities start to degrade very much with increasing Z since the clipping losses augment 
exponentially over the mirror transmissivity losses. We find for instance that using N = 10 
we can afford gaps up to A = 4 m and still keep the minimum g scns (which we take here 
as an indicator of the best attainable sensitivity) comparable to the N = 12 configuration 
with a maximum of A = 2 m gaps. 

In summary, a whole LSW program could be: find the optimal number of magnets with 
the requirements on the best attainable power build-up and magnet aperture. Perform the 
experiment with all possible wiggler configurations. Then, change the gap size to a value 
that does not spoil the whole sensitivity too much and perform experiments in all the 
possible configurations again, this would broaden the sensitivity peaks of the different 
wiggler configurations. The results so far still present the uncomfortable dips caused by 
the single magnet form factor. In order to cope with this, shift all the dips by half an 
oscillation length by introducing a buffer gas as in (51] and perform the experiments above 
again. 

For completeness, we show in Fig. (J7J) a similar optimization analysis for HERA and 
Tevatron magnets. Finally, in Fig. ([S]) we have compared the three types of magnets, in an 
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Figure 6: The quantity l/(L/3 1 / 2 ) to minimize in the design of a LSW experiment looking for ALPs 
(in arbitrary units; and for = 0) as a function of the total length Z of the generation or regeneration 
cavity. We have particularized for the case of the LHC with a = 28 mm. The solid line shows the optimal 
case in which the whole cavity length, Z, can be filled with magnets Z = L. The black dots indicate 
the points in which L is an integer multiple of the LHC magnetic length I = 14.3 m, i.e. the possible 
realistic configurations with negligible gap size. LEFT: The red points show the configurations when the 
gap between the magnets increases to A = 1 m, (Z = N(£ + A)). The dotted line indicates an ideal case in 
which roundtrip losses (including clipping) would be negligible, So — or a = oo. The dashed line assumes 
So = 62- RIGHT: The thin lines emerging from each gaplcss configuration to the right and up represent 
configurations with gap size A = 1, 2, 3, 4, 5 m. 




Z [m] Z [m] 

Figure 7: As before, the solid line shows the optimal case in which the whole cavity length, Z, can be 
filled with magnets Z = L. The black dots indicate the points in which L is an integer multiple of the 
magnet size, I. The red points show the configurations when the gap between the magnets increases to 
A = 1 m, (Z = N(£ + A)). The dotted line indicates an ideal case in which roundtrip losses would be 
negligible, So — or a = 00. The dashed line assumes So = 62- LEFT: HERA setup. RIGHT: Tevatron 
setup. 

ideal gapless situation. As can be inferred from this last plot, the minimum of g sens can be 
shifted nearly by one or two magnets, leading to an erroneous estimation of the sensitivity. 
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Figure 8: Comparison of <7 scn s vs. Z between three possible LSW setups with HERA, LHC, and Tevatron 
magnets, in the ideal scenario of no gap in-between the magnets. 



Sc-dipole magnet 


B(T) 


a (mm) 


L op t (m) 


N 


fens (GeV x ) 


HERA (£ = 8.8 m) 


5.5 


30 


200 


20 


9.8 xlO" 12 


Tevatron (£ = 6 m) 


5 


24 


128 


18 


1.7 xlO" 11 


LHC (£ = 14.3 m) 


9.5 


28 


175 


12 


5.8 xlO -12 



Table 3: Optimal cavity length for three different types of superconducting dipole magnets. The bore 
aperture radius a given assumes straightened magnets. The bounds on g sens have been computed using 
the benchmark values given in Table [T] 

4.2 Example setups 

In Table [3] we show the estimated optimum cavity length for different superconducting 
dipole magnets currently available, the maximum number of magnets that could fix in the 
cavity, N, and the expected sensitivity in coupling constant, based in our analysis from 
the previous section. 

The LHC magnets lead with a coupling constant of 5.8 x 10~ 12 GeV -1 , but it is inter- 
esting that the results obtained for HERA and Tevatron magnets are not that far away. 
Using LHC magnets, and considering a gap between the magnets of 1 m, it would be 
possible to use 12 + 12 magnets, as discussed before, without losing the desired ~ 10 5 . 
For HERA and Tevatron it is possible to arrange optimal configurations of 20 + 20 and 
18 + 18, respectively. In Fig. M we summarize this configurations for the three types of 
superconducting dipole magnets. We have assumed power build-up, power incident, etc, 
as in Tab. [TJ The different "wiggler" configurations have been considered, according to the 
number of magnets from Table |3j 

In the case of massive hidden photons, the inclusion of the gaps does not affect the 
sensitivity of the experiment, since the mixing is not affected by the magnetic field. The 
exclusion plot, considering the optimal cavities for the next generation is shown in FigHTfl 
We have chosen just the HERA configuration for illustration reasons, since all the three 
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Figure 9: Estimate of feasible ALP bounds for the three types of superconducting dipole magnets 
currently available 

configurations are quite similar. 

In the case of massless hidden photon and extra mini-charged particles, the separation 
between the magnets can - in principle - affect the sensitivity of the experiment. In this 
case, the computation of a generic formula for N magnets is more complicated. However, 
we proceeded adding successive gaps to the all order formula, Eq. ([2$. For instance, 
considering just one gap, we have found that their contribution appears in the probability 
of conversion P 7 _>y, at first order in x, as 



With fi 2 = e^AJVj, where AJVj accounts for the complex refractive index induced by the 
creation of a MCP pair and / is a function in the given parameters. 

Therefore, for weak mixing - as it is expected - it is possible to show that increasing 
the number of gaps does not change this kind of suppression. The exclusion plots expected 
for this type of interaction is shown in Fig. [10] (b). 

5 Conclusions 

The next generation of light shining through walls experiments will most likely use arrays 
of superconducting dipole magnets currently available from the large particle colliders such 
as HERA, Tevatron or LHC. The engineering of such arrays might require the existence of 
gaps of zero magnetic field in the generation and regeneration regions of the experiment. We 
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Figure 10: a) Exclusion plot for masslcss hidden photons. The red area corresponds to the current bounds 
from LSW experiments. Blue area corresponds to the expected exclusion plot for the next generation of 
LSW. The probability for the last has been computed according to Eq. (JTTJ) . b) Exclusion plot for hidden 
photons + MCPs. The red area corresponds to the latest bound and the blue area corresponds to the 
expected sensitivity of the next generation of LSW experiments, according to Table [3] We have assumed 
e/, = e. For both plots, we have assumed the optimal cavity of the HERA (200 m) and the benchmark 
input of TableHJ 



have shown that the presence of these gaps affects the photon-H-ALP oscillation probability 
and how to take advantage of them to improve the oscillation probability in certain regions 
of the ALP parameter space where a gapless configuration would be disastrous. This 
technique is similar to introducing phase shift plates or a buffer gas in the oscillation 
regions, but benefits from the absence of optical losses, since no dispersive element is 
introduced in the path of the light. 

In order to maximize the sensitivity those experiments will have to use optical cavities 
encompassing each oscillation region. The finite aperture of the available magnets poses 
constraints on the transversal extent of the optical modes that can be used in the cavities 
and this in turn limits the length of the whole experiment by the beam's divergence. Using 
magnets with the largest bore aperture is therefore desiderable. We have taken all these 
factors into account and proposed optimized setups for the use of HERA, Tevatron or LHC 
magnets, cf. Tab. [3j The reach of these experiments is shown in Figs. [9] and [10] for ALPs 
and other WISPs such as hidden photons and minicharged particles. In the case of ALPs, 
it constitutes a magnificent improvement upon the currently experimental limits pervading 
even the current astrophysical limits and reaching a region of parameter space (g ~ 10~ n 
GeV -1 at small masses) which has been linked to the resolution of a number of recently 
raised astrophysical conundrums. The realization of such an LSW experiment currently 
seems to be the most promising purely laboratory probe capable of testing the ultralight 
ALP hypothesis. Moreover, we believe it can lead the way towards experiments sensitive 
to a ~meV mass QCD axion. 
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